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Abstract. We show two examples of noise-induced synchronization. We study a 1-d 
map and the Lorenz systems, both in the chaotic region. For each system we give 
numerical evidence that the addition of a (common) random noise, of large enough 
intensity, to different trajectories which start from different initial conditions, leads 
eventually to the perfect synchronization of the trajectories. The largest Lyapunov 
exponent becomes negative due to the presence of the noise terms. 

I INTRODUCTION 

The issue of whether chaotic systems can be synchronized by common random 
noise sources has attracted much attention recently [1-6]. It has been reported 
that for some chaotic maps, the introduction of the same (additive) noise in inde- 
pendent copies of the same map could lead (for large enough noise intensity) to a 
collapse onto the same trajectory, independently of the initial condition assigned to 
each of the copies [1]. This synchronization of chaotic systems by the addition of 
random terms is a remarkable and counterintuitive effect of noise. Nowadays, some 
contradictory results exist for the existence of this phenomenon of noise-induced 
synchronization. It is purpose of this paper to give explicit examples in which it is 
shown that one can indeed obtain such a synchronization. Moreover, the examples 
open the possibility to obtain such a synchronization in electronic circuits, hence 
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suggesting that noise-induced synchronization of chaotic circuits can indeed be used 
for encryption purposes. 

Although the issue of which is the effect of noise in chaotic systems was considered 
at the beginning of the 80's [7], to our knowledge the first atempt to synchronize 
two chaotic systems by using the same noise signal was considered by Maritan and 
Banavar [1]. These authors analysed the logistic map in the presence of noise: 



where ^„ is the noise term, considered to be uniformly distributed in a symmetric 
interval [—W, +W]. They showed that, if W was large enough (i.e. for a large noise 
intensity) two different trajectories which started with different initial conditions 
but used otherwise the same sequence of random numbers, would eventually coin- 
cide into the same trajectory. This result was heavily critisizcd by Pikovsky [8] who 
argued that two systems can synchronize only if the largest Lyapunov exponent is 
negative. He then shows that the largest Lyapunov exponent of the logistic map 
in the presence of noise is always positive and concludes that the synchronization 
is, in fact, a numerical effect of lack of precision of the calculation. Furthermore, 
Malescio [2] pointed out that the noise used to simulate Eq.(l) in [1] was not re- 
ally symmetric. This is because the requirement a;„ G (0, 1), Vn, actually leads to 
discard those values for the random number which do not fulfill such condition. 
The average value of the random numbers which have been accepted is different 
from zero, hence producing an effective biased noise, i.e. one which does not have 
zero mean. The introduction of a non-zero mean noise means that we are altering 
essentially the properties of the deterministic map. 

Noise induced synchronization has been since studied for other chaotic systems 
such as the Lorenz model [1,2] and the Chua circuit [3,5]. Synchronization of 
trajectories starting with different initial conditions but using otherwise the same 
sequence of random numbers was observed in the numerical integration of a Lorenz 
system in the presence of a noise distributed uniformly in the interval [0,14^^], 
i.e. again a noise which does not have a mean of zero [1]. Other detailed studies 
[2] also conclude that it is not possible to synchronize trajectories in a Lorenz 
system by adding an unbiased noise. Similarly, the studies of the Chua circuit 
always conclude that a biased noise is needed for synchronization [5]. Therefore a 
widespread belief exists that it is not possible to synchronize two chaotic systems 
by injecting the same noisy signal to both of them. However, in this paper we 
give numerical evidence that it is possible to synchronize two chaotic systems by 
the addition of a common noise which is Gaussian distributed and not biased. We 
analyse specifically a 1-d map and the Lorenz system, both in the chaotic region. 
The necessary criterion introduced in ref. [8] is fully confirmed and some heuristic 
arguments are given about the general validity of our results. Finally, we conclude 
with some open questions relating the general validity of our results. 





II RESULTS 



The first example is that of the map: 



(2) 



where ^„ is a set of uncorrelated Gaussian variables of zero mean and variance 1. 
We use exphcitely 



We plot in Fig. (la) the bifurcation diagram of this map. In the noiseless case, we 
can see the typical windows in which the system behaves chaotically. The associated 
Lyapunov exponent, A, in these regions is positive. For instance, for uj = 0.3 (the 
case we will be considereing throughout the paper) it is A ~ 0.53. In Fig. (lb) we 
observe that the Lyapunov exponent becomes negative for most values of lu for 
large enough noise level e. Again for u — 0.3 and now for e = 0.2 it is A = —0.17. 

For the noiseless case, it is A > and trajectories starting with differential 
initial conditions, obviously, remain different for all the iteration steps, see Fig. (2a). 
However, when moderated levels of noise (e > 0.1) are used, A becomes negative and 
trajectories starting with different initial conditions, but using the same sequence 
of random numbers, synchronize perfectly, see Fig. (2b). 

According to [8], convergence of trajectories to the same one, or loss of memory 
of the initial condition, can be stated as negativity of the Lyapunov exponent. The 
Lyapunov exponent of a map Xn+i = F{xn) is defined as 



It is the average of (the logarithm of the absolute value of) the successive slopes 
F' found by the trajectory. Slopes in [—1, 1] contribute to A with negative values, 
indicating trajectory convergence. Larger or smaller slopes contribute with positive 
values, indicating trajectory divergence. Since the deterministic and noisy maps 
satisfy F' — f one is tempted to conclude that the Lyapunov exponent is not 
modified by the presence of noise. However, there is noise-dependence through the 
trajectory values Xj, i = 1, ...,N. In the absence of noise, A is positive, indicating 
trajectory separation. When synchronization is observed, the Lyapunov exponent 
becomes negative, as required by the argument in [8]. 

By using the definition of the invariant measure on the attractor, or stationary 
probabihty distribution Pst{x), the Lyapunov exponent can be calculated also as 




(3) 




(4) 



A= (log|F'(a;)|) = (log |/'(a;)|) = J Pst{x)log\f{x)\dx 



(5) 



Here we see clearly the two contributions to the Lyapunov exponent: although 
the derivative f'{x) does not change when including noise in the trajectory, the 



stationary probability does change (see Fig. 3), thus producing the observed change 
in the Lyapunov exponents. Synchronization, then, can be a general feature in 
maps which have a large region in which the derivative is smaller than one. 

Noise will be able, then , to explore that region and yield, on the average, a negative 
Lyapunov exponent. 

The second system we have studied is the well known Lorenz model with random 
terms added [9,1]: 

x^p{y- x) 

y — —xz + rx — y + (6) 
z — xy — bz 

^ now is white noise: a Gaussian random process of mean zero and delta correlated, 
(C(^)C(O) = S{t-t'). We have used p = 10, 6 = 8/3 and r = 28 which, in the deter- 
ministic case, e = are known to lead to a chaotic behavior (the largest Lyapunov 
exponent is A 0.9 > 0). We have integrated numerically the above equations 
using the Euler method with a time step At = 0.001. For the deterministic case, 
trajectories starting with different initial conditions arc completely uncorrelated, 
see Fig. (4a). This is also the situations for small values of e. However, when using 
a noise intensity e — AO the noise is strong enough to induce synchronization of 
the trajectories. Again the presence of the noise terms makes the largest Lyapunov 
exponent become negative (for e = 40 it is A ^ —0.2). As in the example of the 
map, after some transient time, two different evolutions which have started in com- 
pletely different initial conditions synchronize towards the same value of the three 
variables (see Fig. (4b) for the z coordinate). One could argue that the intensity 
of the noise is very large. However, the basic structure of the "butterfly" Lorenz 
attractor remains unchanged as shown in Fig. (5). 

In conclusion, we have shown that it is possible for noise to synchronize trajec- 
tories of a system which, deterministically, is chaotic. The novelty of our results 
is that the noise used in the two examples, a 1-d map and the Lorenz system, is 
unbiased, i.e. has always zero mean. 

There still remain many open questions in this field. They involve the develop- 
ment of a general theory, probably based in the invariant measure, that could give 
us a criterion to determine the range of parameters (including noise levels) for which 
the Lyapunov exponent becomes negative, thus allowing synchronization. Another 
important question relates the structural stability of this phenomenon. Any prac- 
tical realization of this system can not produce two identical samples. If one wants 
to use sthocastic synchronization of electronic emitters and receivers (as a means 
of encryption) one should be able to determine which is the allowed discrepancy 
between circuits before the lack of synchronization becomes unacceptable. 
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FIGURE CAPTIONS 

FIGURE 1: (a) Bifurcation diagram of the map given by Eqs.(2) and (3) in 
the absence of noise terms, (b) Lyapunov exponent for the noiseless map (e = 0, 
continuous line) and the map with a noise intensity e = 0.1 (dotted line) and e — 0.2 
(dot-dashed line). 

FIGURE 2: Plot of two realizations x^^\ x^"^^ of of the map given by Eqs. (2) 
and (3). Each realization consists of 10,000 points which have been obtained by 
iteration of the map starting in each case from a different initial condition (100,000 
initial iterations have been discarded and are not shown). In figure (a) there is 
no noise, e = and the trajectories arc independent of each other. In figure (b) 
we have use a level of noise e = 0.2 producing a perfect synchronization (after 
discarding some initial iterations). 

FIGURE 3: Plot of the stationary distribution for the map given by Eqs. (2) 
and (3) in the (a) deterministic case e = 0, and (b) the case with noise along the 
trajectory, e = 0.2. 

FIGURE 4: Same than figure (1) for the z variable of the Lorenz system, Eqs. (6) 
in the (a) deterministic case e = and (b) e = 40. Notice the perfect synchroniza- 
tion in case (b). 

"Butterfly" attractor of the Lorenz system in the cases (a) of no noise e = and 
(b) e = 40. 
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